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We perform spatially dependent tuning of a GalnP photonic crystal cavity using a continuous wave violet 
laser. Local tuning is obtained by laser heating of the photonic crystal membrane. The cavity resonance shift 
is measured for different pump positions and for two ambient gases: He and N 2 . We find that the width 
of the temperature profile induced in the membrane depends strongly on the thermal conductivity of the 
ambient gas. For He gas a narrow spatial width of the temperature profile of 2.8 /xm is predicted and verified 
in experiment. 


Photonic crystal (PhC) cavities are widely studied 
because of their fascinating applications 1-3 . Arrays of 
PhC cavities can form coupled resonator optical waveg¬ 
uides (CROW), which are very promising for slow light 
applications 4 and the study of light localization 5 . Var¬ 
ious fabrication imperfections can cause a disorder in a 
CROW structure which leads to a detuning of cavities 
from the intended resonance and reduces the waveguide 
throughput and bandwidth. Tuning each cavity indepen¬ 
dently can restore cavities in resonance and counteract 
the disorder. There is a variety of methods to change the 
refractive index of a PhC cavity, including free-carrier 
injection 6 , the nonlinear Kerr effect 7 , thermal effects 8,9 , 
oxidation 10-12 and chemical processes 13 . Of these meth¬ 
ods, thermal tuning through local laser heating is easy, 
reversible and can give a steady-state control of the res¬ 
onance properties of the system. However due to the 
diffusion of heat in the PhC membrane thermal control 
of one cavity will affect the neighbor cavities. The width 
of the temperature profile is determined by the sample 
material and the surrounding media. Therefore one can 
expect that by carefully selecting the sample material 
and ambient medium one can control the width of the 
temperature profile. 

In this work we use the semiconductor alloy 
Gao. 51 Ino. 49 P as a sample material and two gases, nitro¬ 
gen and helium, as a surrounding media. Gao. 51 Ino. 49 P 
has a thermal conductivity 14 of 4.9 W/(m-K) which 
is quite small in comparison to other semiconductor 
materials 15 . The thermal conductivity of gases is often 
assumed to be negligible compared to semiconductors, 
such as Si. However, the effect of the gas on the width 
of the thermal profile depends strongly on the ratio of 
the thermal conductivity of the gas and the semiconduc¬ 
tor. Helium has a high thermal conductivity 16 of 0.153 
W/(m-K), which is more than 6 times higher than of 
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nitrogen 16 (0.024 W/(m-K)) and only 32 times smaller 
that that of Gao. 51 Ino. 49 P. Therefore the combination of 
Gao. 51 Ino. 49 P and He should have a high thermal ex¬ 
change efficiency and small width of the temperature 
profile in comparison with other materials. To investi¬ 
gate the width of the temperature profile in PhC mem¬ 
branes we measured the response of the resonance of a 
HO-type cavity to a spatially scanned continuous wave 
(CW) heating laser focused on the membrane. 




FIG. 1. (a) Optical microscope image of a GalnP PhC (red 
line) with an HO cavity (red dotted ellipse) in the center, (b) 
The pump-probe setup. 

The PhC HO nanocavity is implemented in a triangular 
lattice of air holes with a period a=505 nm and a hole 
radius 0.24 a. The hole shift for the HO cavity is 0.16 a. 
A subharmonic hole structure is induced to enhance the 
radiation from the cavity in a vertical direction 17 . The 
cavity is made of a suspended membrane of Gao. 51 Ino. 49 P 
with a thickness of 180 nm. The membrane is covered 
with 30 nm of Si 3 N 4 on the top using a PECVD technique 
to tune the cavity resonance. The detailed description of 
the fabrication method can be found in Ref. 18. A micro¬ 
scope picture of the cavity with the surrounding PhC is 
presented in Fig. la. The cavity is close to the center of 
the membrane and the IR light is evanescently coupled 
through a PhC waveguide. Mode converters at the end 
of the waveguide decrease the Fabry-Perot interference 
in the waveguide and increase the coupling efficiency 19 . 

The setup is shown in Fig. lb. The CW IR laser 
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probe light was coupled to the PhC waveguide using a 
polarization maintaining lensed fiber with NA of 0.55. 
The out-of-plane scattered light was collected using a 0.4 
NA objective and imaged on an IR CCD camera using 
tube lens L2. The total magnification of the system is 
50x. We recorded spectra by sweeping the wavelength 
of the IR laser and taking IR camera snapshots for each 
wavelength. The oxygen concentration was reduced to 
less than 0.025% to suppress oxidation effects 10-12 . The 
chamber was kept at a slight overpressure (less than 2 
mbar above atmospheric pressure). Thermal tuning was 
performed using a 405 nm CW diode laser, which was 
focused onto the surface of the sample into a spot with 
a FWHM of 0.96 //m using the same objective. The sur¬ 
face of the sample with the pump spot was also imaged 
with a visible range camera using tube lens L3 with a 
system magnification of 27x. The position of the pump 
spot was controlled using a spatial light modulator. The 
pump spot was moved along a line through the center 
of the cavity (see Fig. la) and the resonance spectrum 
of the cavity was measured for a sequence of pump po¬ 
sitions. Any effects of slow oxidation or water coverage 
are eliminated by alternating between reference spectra 
(without pump laser) and signal spectra (with pump). 



The two modes have slightly different spatial positions, 
and the reference position is the center of mode Ml. We 
use the wavelength shift of mode Ml to determine the 
spatially dependent cavity response to the pump. Apart 
from the redshift we observed a reduction of the inten¬ 
sity of the out-of-plane scattered light when the pump 
was positioned close to the cavity. We attribute this re¬ 
duction to the fact that the cavity resonance was tuned 
farther away from the PhC band, thereby reducing the 
evanescent coupling to the waveguide. 

The pump power incident on the sample was 110 fi W 
in N 2 atmosphere. The spatial response curve for N 2 at¬ 
mosphere is presented in Fig. 3a. The maximum redshift 
is 4.5 nm and the FWHM spatial width (A x e ff) of the 
tuning curve is 5.1 /im. In He atmosphere a larger pump 
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FIG. 2. Spectra obtained from the out-of-plane scattered light 
for no-pump reference and for pump positions on top of the 
cavity and 5 fi m away from the center of the cavity in N 2 
atmosphere. Color lines are Lorentzian fits. 

Resonance spectra for the reference and for pump posi¬ 
tions of 0 and 5 /mi relative to the center of the cavity in 
N 2 are presented in Fig. 2. The closer the pump spot is 
placed to the cavity, the stronger the redshift is due to the 
better overlap of the temperature profile and the cavity 
mode. We observe two resonance modes: the high- Q low- 
frequency mode (Ml) and the low-Q high-frequency M2. 
Both resonances have a distinct Lorentzian shape and 
quality factors of about 1000 and 300 are estimated for 
modes Ml and M2 respectively. The calculated loaded 
Q-factor of the mode Ml is about 2500, because the cav¬ 
ity is designed to be overcoupled. The resonance M2 cor¬ 
responds to a higher order mode which is predicted from 
calculations. The mode Ml is much brighter than M2. 


FIG. 3. Redshift of the cavity resonance and AT e // for dif¬ 
ferent pump positions with nitrogen (a) and helium (b) as 
an ambient gas. Red curves represent the redshift calculated 
according to the model. From one measurement to another 
we observed an error less than 3% of the shift for N 2 , and less 
than 1% for He. 

power of 220 /rW was used (see Fig. 3b). The maximum 
observed redshift is 6.2 nm and A x e ff is only 3.7 fin 1 . As 
a result of the high thermal conductivity of the He, both 
the redshift (normalized to power) and A x e f / are smaller 
by 30% compared to N 2 . This thermal profile can be used 
to predict thermal crosstalk in structures with multiple 
cavities such as CROWs. Compensation of crosstalk is 
typically possible when the crosstalk is below about 50%, 
i.e. for elements spaced more than Ax e ff/2. A reduction 
of A x e ff therefore directly translates into a higher pos¬ 
sible integration density of independently tuned cavities. 

We have performed numerical calculations to model 
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the resonance wavelength change of the cavity for dif¬ 
ferent pump positions taking into account material and 
pump properties. For small perturbations of the refrac¬ 
tive index the wavelength change can be expressed as 20 


aa = A^ 4t . 

(n) dJ 


eff 


with 


AT eff = 


J AT(r — r 0 )e(r)|E(r)| 2 rf 3 r 
f e(r)|E(r)| 2 d 3 r 


(1) 

( 2 ) 


where AT e //is temperature averaged over the mode pro¬ 
file e(r)|E(r )| 2 of the cavity, e(r) is the dielectric constant 
of the membrane material, AT = T — To is the local tem¬ 
perature change in the PhC averaged over the membrane 
thickness, ^ is the linear coefficient of the refractive in¬ 
dex temperature response, (n) is the averaged refractive 
index of the membrane material and ro is the position 
of the pump. To = 293.15IF is the ambient temperature. 
The mode profile was calculated using a finite difference 
time domain method 21 . Temperature profiles were cal¬ 
culated in COMSOL using a steady state heat diffusion 
model with a constant source term assuming that ther¬ 
mal conductivities are not temperature dependent within 
the relevant temperature range: 


K(x,y,z)V 2 T(x,y,z) + W(x,y,z) = 0, (3) 


where k(x, y , z) is the thermal conductivity and W is a 
source term inside the membrane written as a sum of two 
terms W = We + Wr, due to the excess energy of the 
carriers and their recombination respectively. 

We observe photoluminescence (PL) of Gao.51Ino.49P 
at 670 nm caused by the pump light. The bandgap of 
Gao.51Ino.49P is 1.85 eV and the energy of a 405 nm 
photon is 3.06 eV. The excess energy term (We) is re¬ 
leased within about a picosecond upon optical excitation 
and therefore has the spatial profile of the pump spot 
(a Gaussian width a E = 0.41 /mi) and corresponds to a 
fraction f E = 0.4 of the initial power. The recombina¬ 
tion term (Wr) has the profile of the PL spot, a Gaussian 
width <jr = 0.50 /im, and it takes a fraction /r = 0.6 of 
the photon energy. We assume that all absorbed light is 
converted into heat, because the estimated radiative car¬ 
rier recombination efficiency is less than 1 %. The excess 
energy term is expressed as 

W E (x,y,z) = f E P{ ^J i] ae~^ + y 2 ^e- az . (4) 

Here R is the reflection coefficient of the membrane es¬ 
timated to be 9.5%, a is the absorption coefficient of 
Gao.51Ino.49P, and P is the total power of the pump 
beam. The recombination heating term has the same 
Gaussian shape with the appropriate fraction and width 
inserted. The size of the air gap between the membrane 
and GaAs substrate for our samples is equal to 1.5 /im, 
the size of the air layer above the membrane was taken 
to be 10 fim. All external boundaries are Dirichlet type 
with temperature Tq. The photonic crystal is assumed 


to be a rectangle of 40x20 /mi. Dirichlet boundaries 
are considered because the membrane is surrounded from 
four sides by the bulk material, which acts like a perfect 
heat sink. Thermal conductivities for gases were taken 
from Ref. 16. The coating layer made of SiaN 4 has a 
thermal conductivity 22 of 24.5 W/(m-K). The thermal 
conductivity of the membrane was taken as volume av¬ 
eraged between holes of the PhC lattice and two layers 
of SiaN 4 and Gao. 51 Ino. 49 P. The absorption coefficient 23 
of Gao. 51 Ino. 49 P at 405 nm is 24.5 (/im ) -1 and the ab¬ 
sorption coefficient 24 of SiaN 4 is less than 0.006 (/mi) . 

As a result nearly all pump light is absorbed by the 
Gao. 51 Ino. 49 P membrane. The temperature profile can 
be calculated without any free parameters. The only pa¬ 
rameter which is unknown for the redshift calculation is 
for Gao. 51 Ino. 49 P- 

The resulting effective temperature profiles are pre¬ 
sented in Fig. 3a,b, and are in excellent agreement with 
the experiment for both nitrogen and helium. The max¬ 
imum temperature rise in case of nitrogen atmosphere is 
29 K, and in case of He it is 37 K, so the ratio of these two 
values normalized to the incident power is equal to 0.63. 
The experimental ratio of redshifts is 0.69. From Fig. 3 
we also estimate the so far unknown ^ for Gao.51Ino.4gP. 
We obtain the value of 3.1 ± 0.5 x 10 _ 4 K _1 . 
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FIG. 4. X-axis crossection of the calculated temperature pro¬ 
file (AT) for different ambient media. The pump focus is 
placed in the center of the cavity. The power of the pump 
light is 220 /rW. The inset shows the 2D averaged tempera¬ 
ture profile for helium. The isotherms have a step of 6.3 K. 

Instead of helium one could use media with a higher 
thermal conductivity. In Fig. 4 calculated temperature 
profiles along the central line of the PhC membrane (see 
Fig. 1) are presented for different surrounding media. 
These are bare temperature profiles, not convolved with 
the cavity mode. The widest temperature profile is found 
in case of vacuum. Water has a thermal conductivity 4 
times higher than He, which results in a significantly nar¬ 
rower temperature profile. We note that low-refractive- 
index liquids such as acetone and isopropanol have ther¬ 
mal conductivities comparable to that of He. Using liq- 
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uids for heat exchange completely changes the optical 
properties of the sample, therefore using a high-thermal- 
conductivity gas is a more appropriate option in many 
cases. 

The general tendency of the width A xth and the peak 
value (AT m ax) of the bare temperature profile is pre¬ 
sented in Fig. 5. In case of low ambient thermal conduc- 



FIG. 5. Calculated Axth (a) and ATmax (b) vs thermal con¬ 
ductivity ratio of surrounding media and Gao. 51 Ino. 49 P. Some 
well-known materials are marked. The power of the pump 
light is 220 /iW. 

tivity heat tends to diffuse inside the membrane, there¬ 
fore A x t h and the peak temperature reach constant val¬ 
ues of 5 /im and 79 K, respectively. For large thermal 
conductivities of the surrounding media all heat diffuses 
rapidly into the ambient media reducing A xth and the 
peak temperature. 

Gas cooling is most relevant for thermal tuning of reso¬ 
nances in membranes made of materials with small ther¬ 
mal conductivities such as Gao.51Ino.49P. In highly con¬ 
ducting materials like Si the thermal profile is signifi¬ 
cantly larger 9 and is typically determined by the material 
boundaries rather than by heat exchange with the gas. 

For time-dependent tuning applications such as a 
switch or a tracking filter the time response is of im¬ 
portance. From time-dependent simulations we estimate 
that the 10% to 90% temperature rise time in a nitrogen 
atmosphere is about 6 [is, but in helium it reduces to 
2 [is. Therefore, helium not only helps to improve the 
spatial resolution of thermal tuning but also one makes 
the tuning response about three times faster. Similarly, 
it should reduce thermal optical nonlinearities that result 
from absorption of cavity-enhanced probe light. 

In conclusion, we measured the spectral tuning curve 


of a thermally tuned cavity for different surrounding me¬ 
dia and found an excellent agreement between experi¬ 
ment and our numerical calculations. Our results and 
our calculations are of importance for applications such 
as thermal control of CROWs where the width of the 
temperature profile should be minimized. Also our cal¬ 
culations offer insight in the role of ambient materials for 
any applications that involve local laser heating. 
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